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We present results related to determination of the Unitarity Triangle angle cj>3. 



1 Introduction 

Determinations of the Cabibbo-Kobayashi-Maskawa (CKM) 1 matrix elements provide important 
checks on the consistency of the Standard Model and ways to search for new physics. Various 
methods using CP violation in B — > DK decays have been proposed 2 ' 3 > 4 ' 5 > 6 to measure the 
Unitarity Triangle angle fa . These methods are based on two key observations: neutral D° 
and D° mesons can decay to a common final state, and the decay B + — > D^K + can produce 
neutral D mesons of both flavors via b — ► cus (Fig. la) and b — > ucs (Fig. lb) transitions, with 
a relative phase 6+ between the two interfering amplitudes that is the sum, 5 + 03, of strong 
and weak interaction phases. For the charge conjugate mode, the relative phase is 9- = 5 — fa. 

The results are based on a data sample containing 275 million BB pairs, collected with the 
Belle detector 7 at the KEKB asymmetric energy e + e~ collider operating at the T(45) resonance. 
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Figure 1: Feynman diagrams of (a) dominant B + — > D°K + and (b) suppressed B + — > D°K + decays 



Table 1: Yields and asymmetries obtained for Dh modes. 



Dfir 19283 ± 150 9690 ± 104 9521 ± 103 -0.01 ± 0.01 

Di?r 2183 ± 55 1051 ± 36 1132 ± 37 0.04 ± 0.02 

P- 2 vr 2413 ± 93 1178 ± 43 1226 ± 43 0.02 ± 0.03 

D f K 1031 ± 39 484 ± 27 549 ± 28 0.06 ± 0.04 

D\K 114 ±21 49 ± 15 63 ± 14 0.07 ± 0.14 

D 2 K 167 ±21 94 ± 17 75 ± 16 -0.11 ± 0.14 



2 B ± -» DcpK* and D* CP K^ 

Recent theoretical studies on 5 meson dynamics have demonstrated a method to access 4>z using 
the process B~ -> dW°^W _ 2 . When a D° is reconstructed CP eigenstate, the 6 — > c 
and b —> u processes interfere. This interference leads to direct CP violation. To extract </> 3 and 
assuming no D°-D° mixing, some necessary observables sensitive to CP violation are: 



Al,2 = 



B(B~ -> D 1>2 K~) - B(B+ -» D 1>2 K + ) 2r sin 8' sin </> 3 



B(B- -> D 1)2 K~) + B{B+ -» D^ 2 K+) 1 + r 2 + 2r cos <$' cos </> 3 



7^ 2 = -^ = l + r ±2r cos 5 cos 03, 5 = j tf + ff for ^ 
where the ratios P Dl > 2 and P D ° are defined as 



R D ° 



B{B~ -> g^^-J + ^ -Pi,2^ + ) 

6(5- -> Di, 2 7r-) + B(B+ -> L> 1;2 7r+) 
-» D°K~) + B{B + -» Z)°^+) 



where L>i and L> 2 are CP-even and CP-odd eigenstates respectively. The asymmetries *4i and 
A 2 have opposite signs. The ratio r is defined as r = \A(B~ — > D°K~)/A(B~ — ► and 
is the ratio of the two tree diagrams shown in Fig. 1 where 5 is their strong-phase difference. 
The size of the ratio r governs the magnitude of the maximum possible CP asymmetry; this 
ratio is suppressed by both CKM (~ 0.45) and color (~ 0.40) factors. The asymmetries and 
double ratios can be calculated for D* in a similar manner. The analysis is described in detail 
elsewhere 8 . 

Fig. 2 shows the AE distributions for B ± — > DK ± events. Table 1 summarizes the yields 
from AE fit and the corresponding asymmetries with statistical errors. In the control samples, 
no large deviation from is seen. The modes of interest are D\K and D 2 K where the B + and 
B~ events are used to calculate asymmetries and double ratios. 
The final asymmetries for A\ and A 2 are found do be 

Ai = 0.07 ± 0.14 (stat) ± 0.06 (sys) 
A 2 = -0.11 ± 0.14 (stat) ± 0.05 (sys) 

agreeing with theoretical expectations where they should have opposite signs. The double ratios: 

Pi = 0.98 ± 0.18 (stat) ± 0.10 (sys) 
P 2 = 1.29 ± 0.16 (stat) ± 0.08 (sys) 
Fig. 3 shows the AE distributions for P+ — > D*K ± events. Table 2 contains the yields of 
the distributions with statistical errors and asymmetries. The statistical significance of D\K 
and D 2 K signals are 5.6 and 4.5 respectively. 
Asymmetries were found to be: 
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Figure 2: AE distributions for (top left) B + -> , (top right) B~ 

(bottom right) B~ -> D 2 A' _ . 



DiA" - , (bottom left) B+ -> 




Figure 3: AE distributions for (left) B ± -> DfK*, (right) B ± -> D^ ) K ± . 



Table 2: Yields and asymmetries obtained for D*h modes. 

D*tt 5762 ± 101 2681 ± 74 2594 ± 74 -0.02 ± 0.02 

Dfir 795 ± 41 399 ± 24 397 ± 23 0.00 ± 0.04 

D|tt 715 ± 37 415 ± 35 301 ± 33 -0.16 ± 0.07 

D*K 284 ± 23 158 ± 16 127 ± 16 -0.11 ± 0.08 

D\K 56 ± 11 33 ± 8 19 ± 8 -0.27 ± 0.25 

D* 2 K 33 ± 10 13 ± 6 22 ± 7 0.26 ± 0.26 



Ai = -0.27 ± 0.25 (stat) ± 0.04 (sys) 

A 2 = 0.26 ± 0.26 (stat) ± 0.03 (sys) 
where the systematic errors were calculated in a similar way to the Dh case. Double ratios 
found are: 

Ri = 1.43 ± 0.28 (stat) ± 0.06 (sys) 
R 2 = 0.94 ± 0.28 (stat) ± 0.06 (sys) 

In summary, the partial rate asymmetries A12 are measured for the decays — > D^jpK^ 
and are consistent with zero. A first observation is seen for D\K and D\K. 

3 Measurement of <fi 3 with Dalitz Plot Analysis of B^ 1 — »■ D^K^ 1 Decay 

Recently, three body final states common to D° and D°, such as Kstt + tt~ 6 , were suggested 
as promising modes for the extraction of fo. This method is based on two key observations: 
neutral D° and D° mesons can decay to a common final state such as K s ir + ir~ , and the decay 
can produce neutral D mesons of both flavors via b — > cus and 6 — > ucs 
transitions, where the relative phase 0+ between the two interfering amplitudes is the sum, 
5 + (^3, of strong and weak interaction phases. In the charge conjugate mode, the relative phase 

= 5 — 4>3, so both phases can be extracted from the measurements of such B decays and their 
charge conjugate modes. The phase measurement is based on the analysis of Dalitz distribution 
of the three body final state of the D meson. The analysis is described in detail elsewhere 9 . 

The Dalitz plots of D decaying to K s it + -it~, which contain information about CP violation 
in B decays, are fitted for B~ and B + data sets. A combined unbinned maximum likelihood 
fit to the B + and B~ samples with r, <j>3 and 5 as free parameters yields the following values: 
r = 0.25 ± 0.07, 3 = 64° ± 15°, S = 157° ± 16° for the B ± -» DK ± sample and r = 0.25 ± 0.12, 
(f> 3 = 75° ±25°, 5 = 321° ±25° for the B ± -» b*K± sample. The errors quoted here are obtained 
from the likelihood fit. These errors are a good representation of the statistical uncertainties 
for a Gaussian likelihood distribution, however in our case the distributions are highly non- 
Gaussian. In addition, the errors for the strong and weak phases depend on the values of the 
amplitude ratio r (e.g. for r = there is no sensitivity to the phases). A more reliable estimate 
of the statistical uncertainties is obtained using a large number of MC pseudo-experiments as 
discussed below. 

We use a frequentist technique to evaluate the statistical significance of the measurements. 
To obtain the probability density function (PDF) of the fitted parameters as a function of 
the true parameters, which is needed for this method, we employ a "toy" MC technique that 
uses a simplified MC simulation of the experiment which incorporates the same efficiencies, 
resolution and backgrounds as used in the data fit. This MC is used to generate several hundred 
experiments for a given set of r, 9 + and 9- values. For each simulated experiment, Dalitz plot 
distributions are generated with equal numbers of events as in the data, 137 and 139 events for 
B~ and B + decays, correspondingly, for B^ — > DK^ mode and 34 and 35 events for B~ and 
B + for B^ — > D*K ± mode. The simulated Dalitz plot distributions are subjected to the same 
fitting procedure that is applied to the data. This is repeated for different values of r, producing 
distributions of the fitted parameters that are used to produce a functional form of the PDFs of 
the reconstructed values for any set of input parameters. 

The confidence regions for the pairs of parameters (<ps,5) and ((f>3,r) are shown in Fig. 4 
(B^ — > DK^ mode) and Fig. 5 (B^ — > D*K^ mode). They are the projections of the corre- 
sponding confidence regions in the three-dimensional parameter space. We show the 20%, 74% 
and 97% confidence level regions, which correspond to one, two, and three standard deviations 
for a three-dimensional Gaussian distribution. 

For the final results, we use the central values that are obtained by maximizing the PDF and 
the statistical errors corresponding to the 20% confidence region (one standard deviation). Of 
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Figure 4: Confidence regions for the pairs of parameters (a) (r, (pa) and (b) (<f>3,5) for the B ± — ► _DA' ± sample. 
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Figure 5: Confidence regions for the pairs of parameters (a) (r, 4>z) and (b) (03,(5) for the B ± — > D*K ± sample. 

the two possible solutions (03, 8 and 03 + 180°, 8 + 180°) we choose the one with < 03 < 180°. 
The final results are 

r = 0.21 ± 0.08 ± 0.03 ± 0.04, 3 = 64° ± 19° ± 13° ± 11°, 8 = 157° ± 19° ± 11° ± 21° (1) 

for the i? ± — ► DK^ mode and 

r = 0.12ig;Jf ± 0.02 ± 0.04, 3 = 75° ± 57° ± 11° ±11°, 5 = 321° ± 57° ± 11° ± 21° (2) 

for the B ± — > D*K ± mode. The first, second, and third errors are statistical, systematic, and 
model dependent errors. 

The significance of CP violation is 94% for the — > DK^ sample and 38% for — > 

The two events samples, B ± — > Di^ 1 * 1 and B ± — > D*K ± , are combined in order to obtain a 
more accurate measurement of 03. The 03 result from the combined analysis is 

3 = 68° t\fo ± 13° ± 11°, (3) 

where the first error is statistical, the second is experimental systematics, and the third is 
model uncertainty. The two standard deviation interval including the systematic and model 



Table 3: Summary of the results. For the B~ — » D sup K~ signal yield, the peaking background contribution 
has been subtracted. The first two errors on the measured production branching fractions are statistical and 
systematic, respectively, and the third is due to the uncertainty in the B~ — > D{ av h~ product branching fraction 
used for normalization. 



Mode 




Signal Yield 


Statistical 
significance 


Measured product 
branching fraction 


Upper limit 
(90%C.L.) 


B~ -> 
B~ -> 
B~ 
B~ -> 


D snp K- 

D mp TT~ 

D iav K- 


8.5^ 
28.5 _7 4 
376.0 

8181.9 1 9 ^ 


2.3cr 
6.4a 


(3.2 +^±0.2 ±0.5) x 10" 7 
(6.6 ±i;? ±0.4 ±0.3) x 10~ 7 


6.3 x 10~ y 




uncertainties is 22° < </>3 < 113°. The statistical significance of CP violation for the combined 
measurement is 98%. 

4 Study of the Suppressed Decays B — * [K + ir ~]dK~ and B~ — > [K + Tt ~]z)7r _ 

As noted by Atwood, Dunietz and Soni (ADS) 5 , CP violation effects are enhanced if the final 
state is chosen so that the interfering amplitudes have comparable magnitudes; the archetype 
uses B~ — > [K + 7t~]dK~ , where [-fT + 7r~]£> indicates that the K + tt~ pair originates from a 
neutral D meson. The analysis is described in detail elsewhere 10 . 
The ratio of branching fractions is defined as 

= B(B- -» D sup h-) _ N Dsnph - /e Dsuph - 
Dh B(B~ - Dtohr) N D ^ h -/e D{avh - ' 1 > 

where Nf) sup h (ND fav h) an d eD Bup h- ( e r> fav fc-) are * ne number of signal events and the reconstruc- 
tion efficiency for the decay B — > D sup h~ (B~ — > Df av /i~), and are given in Table 3. 
The ratios Roh are calculated to be 

R DK = (2.3i^j(stat) ±0.1(syst)) x 10~ 2 , 
R Dn = (3.5l^(stat) ±0.2(syst)) x 10~ 3 . 




Table 4: Signal yields and partial rate asymmetries. 



Mode 


N(B~) 


iV(S+) 




B -» D sup K 

B -> AupTT 


o 9+5.U 

18.8 i|| 


in 1 + 5 - 5 
iu - i -4.8 


0.88 +^±0.06 
0.30 toil ±0.06 



Since the signal for B — ► L> sup .K" is not significant, we set an upper limit at the 90% confidence 
level (C.L.) of R DK < 4.4 x 10~ 2 . 

The product branching fractions for B~ 



B{B~ 



D sup h' 



)=B(B- 



are determined as 
Dfavh~) x R Dh , 



(5) 



and are given in Table 3. A third uncertainty arises due to the error in the branching fraction 
of B~ — ► Df av h~, which is taken from 11 . The uncertainties are statistics-dominated. For the 
B~ — > D sup K~ branching fraction, we set an upper limit at the 90% C.L. of B{B~ — > D sup K~) < 
6.3 x 10~ 7 . For B — > Z? sup 7r~, our measured branching fraction is consistent with expectation 
neglecting the contribution from B~~ — ► D°tt~ . 
The ratio Rdk is related to 03 by 



-Rdx = r 2 B +r 2 D + 2r B r D cos 3 cos 5, 



where 11 



r D = 











^D°K-) 




A(£>° - 


-» 2TV") 




A(D° - 







0.060 ± 0.003, 



and <5b and <5d are the strong phase differences between the two 1? and D decay amplitudes, 
respectively. Using the above result, we obtain a limit on tb- The least restrictive limit is 
obtained allowing ±lcr variation on rjj and assuming maximal interference (03 = 0°,<5 = 180° 
or 03 = 180°, 5 = 0°) and is found to be r B < 0.27. 



We search for partial rate asymmetries Adh in B^ — > D snv ,h^ decay, fitting the B and B + 
yields separately for each mode, where Adh is determined as 

A = B{B- -> D sup h~) - B(B+ -> ZJ sup /i+) 

The peaking background for B T — > Ampi^ is subtracted assuming no CP asymmetry. The fit 
results are shown in Fig. 7 and Table 4. We find 

_4 DK = 0.88l^(stat) ±0.06(syst), 
Adit = 0.30±g;||(stat) ±0.06(syst). 

In summary, we observe B~ — > D sup 7r~ for the first time, with a significance of 6.4<j. The 
size of the signal is consistent with expectation based on measured branching fractions 11 . The 
significance for B~ — > D sup K~ is 2.3<r and we set an upper limit on the ratio of B decay 
amplitudes re < 0.27 at 90% confidence level. 
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